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SUMMARY 

An  introduction  is  given  to  the  application  of  boundary  Integral 
equation  methods  to  the  areas  of  wave  propagation,  radiation  and 
scattering.  Such  applications  were  among  the  first  for  this  method 
which  made  use  of  the  power  of  modern  digital  computers;  however, 
there  is  a long  history  of  integral  methods  in  this  field.  The 
emphasis  here  is  on  fundamental  concepts  coupled  with  illustrations 
of  the  application  of  this  method.  The  breadth  of  the  applications 
of  this  approach  precludes  any  attempt  at  overall  compl-.^.is  and 
many  details  are  left  to  an  extensive,  but  undoubtedly  incomplete, 
list  of  references.  There  is  however  sufficient  material  to  formulate 
the  solution  of  most  problems,  with  some  of  the  inevitable  pitfalls 
indicated. 
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INTRODUCTION 


In  view  of  the  role  of  integral  equations  in  the  development  of 
wave  theory,  it  is  not  strange  that  some  of  the  first  applications  of 
boundary  integral  equation  methods  were  to  examples  in  wave  propagation 
and  scattering.  The  use  of  Green's  functions  in  all  branches  of 
mathematical  physics  is  pervasive  enough  that  individual  references 
are  scarcely  required;  one  readily  available  reference  which  includes 
a wide  range  of  applications  of  Green's  functions  is  Morse  and  Feshbach 
(1953).  A review  of  Green's  function  techniques  used  by  the  mathemati- 
cal physics  community  at  that  time,  as  exemplified  by  this  text,  would 
however  reveal  that  this  use  was  primarily  either  to  obtain  analytical 
(exact)  solutions  to  a fairly  restricted  class  of  problems  (e.g., 
those  simple  geometries  and  boundary  conditions  for  which  a particular 
Green '8  function  could  be  constructed  so  as  to  reduce  Green's  Theorem 
to  a simple  guadrature)  or  to  obtain  asymptotic  expansion  solutions. 

A third  form  of  application  lay  in  the  de  ^lopment  cl  various  integral 
identities,  existence  and  uniqueness  proofs  etc.,  -l.Dse  have  played 
a major  role  in  the  development  of  general  theory. 

The  use  of  boundary  integral  equations  to  solve  wave  problems 
numerically  is  a fairly  recent  development  with  the  first  papers  using 
computer  solutions  based  on  the  boundary  integral  equation  method  for 
transient  acoustic  wave  scattering  appearing  in  the  early  1960 's, 
e.g.,  Friedman  and  Shaw  (1962),  Shaw  and  Friedman  (1962),  followed  by 
a similar  treatment  of  time  harmonic  elastic  and  acoustic  wave  scat- 
tering, Banaugh  and  Goldsmith  (1963).  Prior  to  this  time,  there  was 
of  course  a great  deal  of  work  done  on  these  forms  of  equations,  e.g., 
Rupradse  (1935) , but  these  papers  appear  to  mark  the  shift  to  computer 


solutions.  At  this  same  time,  Hess  (1962)  was  developing  computer 
solutions  to  potential  flow  problems  and  Jaswon  (1963) , Symm  (1963) 
for  other  potential  problems.  Since  that  time,  the  BIB  method  has  been 
employed  (in  some  cases,  rediscovered)  as  a numerical  solution  method 
by  researchers  in  many  areas  of  wave  research.  Before  going  into 
details  of  individual  cases,  it  is  useful  to  separate  these  problems 
into  sub-groups  by  recognising  some  distinctive  characteristics.  One 
separation  ie  that  of  transient  versus  time  harmonic  behavior;  both 
are  of  great  physical  Importance  in  many  areas  of  physics  and  engineer- 
ing but  require  quite  different  formulations  and  solution  methods. 
Another  division  is  that  between  exterior  and  interior  problems;  for 
the  former,  solutions  are  sought  in  an  infinite  domain  exterior  to 
some  prescribed  surface  while  the  latter  requires  solution  in  a bounded 
region  interior  to  such  a surface.  (Ihi6  distinction  blurs  somewhat 
for  surfaces  of  infinite  extent) . Appropriate  boundary  conditions  are 
specified  on  the  surface  in  either  case  with  the  exterior  problem 
requiring  some  further  specif icatiou  of  'behavior  at  infinity*.  A 
further  grouping  may  be  made  on  the  basis  of  the  inhomogeneity  or 


forcing  function.  A problem  driven  by  some  prescribed  boundary  condi- 
tion, e.g.,  a given  displacement  or  velocity  field  on  a specified 
surface  would  be  considered  a mechanical  radiation  problem.  A prob- 
lem wherein  an  obstacle  with  a prescribed  boundary  condition  (usually 
homogeneous)  Interacts  with  some  Incident  wave  field  generated  by 
sources  elsewhere  (frequently  considered  as  plane  waves  coming  from 
'infinity')  would  be  called  a scattering  problem.  Of  Interest  as 
well  Is  the  unforced  or  eigenvalue  problem.  Still  further  distinction 
may  be  made  between  two  and  three  spatial  dimenalonal  problems;  this 


3 


distinction  being  straightforward  in  time  harmonic  problems  but  pre- 
senting a serious  difficulty  in  transient  problems.  Other  categoriza- 
tion may  be  made  on  the  basis  of  the  type  of  boundary  conditions 
prescribed  (Dirichlet,  Neumann,  Robin,  etc.).  Finally, and  possibly 
most  Important  physically,  is  the  particular  medium  through  which  the 
waves  propagate,  e.g.,  acoustic,  elastic,  electromagnetic,  water 
(especially  free  surface  gravity  waves) , and  thus  the  physical  signi- 
ficance of  the  dependent  variables  and  their  boundary  and  initial  con- 
ditions. Fortunately,  the  boundary  integral  equation  approach  is 
applicable  to  all  of  these  'variations  on  a theme' . Its  primary  limi- 
tations are  the  requirement  of  a known  fundamental  solution  for  the 
particular  governing  equation  and  linearity  of  the  equation,  but  not 
necessarily  of  the  boundary  conditions.  Emphasis  here  will  be  given 
to  mechanical  rather  than  electromagnetic  wave  problems;  integral 
equation  approaches  to  the  latter  are  discussed  at  length  by  Mittra 
(1973). 

The  bulk  of  the  research  based  on  boundary  integral  equation  metl.oua 
applied  to  wave  problems  has  been  directed  towards  numerical  solution, 
using  large  scale  digital  computers  and  this  will  form  the  primary 
theme  in  this  chapter.  Since  efficient  algorithms  exist  for  the  solu- 
tion of  systems  of  algebraic  equations,  the  most  common  procedure  is 
to  reduce  the  governing  boundary  integral  equation  to  such  a system  by 
approximating  the  boundary  values  of  the  dependent  variable  by  suitable 
trial  functions  over  selected  boundary  elements.  The  nature  of  the 
boundary  Integral  equation  approach  is  to  reduce  the  number  of  approxi- 
mating dependent  variables,  e.g.,  trial  function  coefficients,  by  only 
requiring  boundary  values  of  the  original  dependent  variable  as  opposed 
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to  other  solution  schemes  which  require  s solution  everywhere  In  the 
domain  of  the  problem.  Such  a procedure  is  particularly  useful  in 
exterior  problems  where  other  numerical  solution  methods  require  some 
'finite'  closing  surface  at  which  their  approximating  mesh  may  be  ter- 
minated, e.g.,  see  Vastano  and  Reid  (1967)  and  Shaw  (1974)  for  a 
direct  comparison  of  exterior  solutions  by  finite  differences  and 
BIE  respectively.  The  price  paid  for  this  reduction  is  a full  coeffi- 
cient matrix  as  opposed  to  the  sparse  matrices  found  in  corresponding 
finite  difference  or  finite  element  methods.  Some  overall  reviews 
of  this  approach  are  available;  Schenck  (1967)  for  time  harmonic 
acoustic  scattering  and  radiation,  Shaw  (1970)  for  transient  and  time 
harmonic  acoustic  scattering  and  radiation,  Chertock  (1971)  and 
Klelnman  and  Roach  (1974)  again  for  Helmholtz  problems  and  Nlttra  (1973) 
for  the  electromagnetic  case.  In  addition,  the  books  by  Kupradze  (1935), 
(1963)  discuss  integral  equation  formulations  for  elastodynamlcs , 
while  Pao  and  Mow  (1973)  give  a section  on  numerical  solution  for  elas- 
tic scattering  problems  using  boundary  integral  equations.  It  is  clear, 
however,  the  boundary  integral  equation  method  is  more  than  a numerical 
technique;  it  is  a reformulation  completely  equivalent  to  any  original 
partial  differential  equation  representation.  As  such,  it  lands  itself 
to  other  forms  of  solution,  l.e.,  analytical  and  approximate  as  wall 
as  numerical.  Analytical  solutions  are  widely  covered  in  classical  mathe- 
matical physics,  e.g.,  Morse  and  Feshbach  (1953),  and  will  not  be 
pursued  further  here.  The  distinction  between  approximate  and  maarlcal 
solutions  is  one  of  semantics  since  clearly  numerical  solutions  are 
also  approximate  in  some  sense.  What  is  meant  here  is  that  approximate 
solutions  will  be  taken  as  global  or  functional  solutions  to  a system 
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of  equations  which  represent  some  form  of  expansion  of  the  original 
equation  in  terms  of  a physical  parameter  which  characterizes  the 
problem.  Numerical  solutions  on  the  other  hand  will  be  solutions  to 
the  original  equation  with  some  piecewise  and  local  approximate  rep- 
resentation of  the  dependent  variable.  Thus  in  the  former,  it  is 
the  governing  equation  which  is  'approximated'  while  in  the  latter  it 
is  the  dependent  variable.  Typical  of  such  approximate  techniques  are 
the  high  frequency  asymptotic  expansions  of  classical  optics,  e.g., 
Northover  (1971).  Also  of  Interest  are  expansions  in  terms  of  geometri- 
cal parameters,  for  example  the  offset  distance  in  shells  of  slightly 
nonconcentric  circular  or  other  simple  boundaries,  Shaw  and  Tai  (1974), 
Shaw  (1978a)  or  the  plate  thickness  in  the  theory  of  thin  elastic 
plates,  Shaw  (1978b).  Examples  of  the  latter  approach  will  be  dis- 
cussed in  a later  section.  Solution  schemes  based  on  Galerkin  or 
related  moment  methods  can  be  employed  as  well  in  the  solution  of 
integral  equations,  e.g..  Delves  and  Walsh  (1974),  Fr~  (1962)  for 
general  theory  and  Sharma  (1967)  for  a specific  example,  indeed  the 
piecewise  numerical  approximation  mentioned  above  could  be  thought  of 
as  a snecial  case  of  a Galerkin  approach  with  particularly  simple 
trial  functions. 

FORMULATION 

Wave  propagation  problems  arise  in  so  many  areas  of  mathematical 
physics  and  engineering  sciences  that  it  is  important  to  concentrate 
on  their  similarities,  which  are  primarily  mathematical,  rather  than 
their  differences,  which  are  primarily  physical.  In  some  cases,  the 
physical  differences  require  a difference  in  formulation  and  these 
will  be  treated  as  they  arise.  In  general  however,  wave  propegetlon 
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probleats  may  be  classified  as  hyperbolic  (elate  dependent)  or  elliptic 
(time  harmonic)  differential  systems . The  most  common  governing 
equations  of  these  forms  are  the  wave  equation 


(1) 


V2  (?,t)  = /c 


and  the  Helmholtz  equation  (with  a time  harmonic  behavior  required) 
i.e.,  <|»(r,t)  - $(r)exp[-iutt] 

(2)  y ! <$>  ( r ) + <}>  ' = ° 

respectively  with  c the  phase  speed  and  k - u/c  the  wave  number. 

The  integral  equation  formulation  for  the  wave  equation,  developed 
in  1882  by  Kirchhoff,  has  been  widely  interpreted  as  a mathematical 
formulation  of  an  even  earlier  physical  statement,  Huygens  principle, 
e.g.,  see  Baker  and  Copson  (1939).  While  the  form  of  Kirchhoff 's 
Integral  formulation  can  be  derived  directly  by  an  application  of 
Green's  Theorem  in  four  dimensional  space  time,  e.g.,  Friedman  and 
Shaw  (1962)  , it  follows  most  readily  albeit  in  a more  restricted  sense 
from  Helmholtz's  integral  formulation;  this  shall  be  the  approach 
here. 

Consider  the  Helmholtz  equation  with  an  inhomogeneous  term, 
applied  to  a domain  D with  a boundary  B upon  which  some  appropriate 
boundary  conditions  on  + (e.g.,  Dlrlchlet , Neumann,  Robin)  are 

specified.  The  fundamental  or  point  source  or  infinite  space  Green's 
function  solution  (these  terms  shall  be  used  here  Interchangeably) 
satisfies  (using  V*  to  represent  the  Lap  laden  in  r^  coordinates) 
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(j)  V„lG(v:;  i'J  + IfG  (■*;*•)  = " +Tr  S(r'r°^ 

and  may  be  found  as 

(«>  &{?;?.)-- G(\?-^)--  «*pCi,tR3/R 

where  R = |r  - rQ|.  The  sign  in  the  exponential  term  is  chosen  to 
insure  outgoing  waves  from  the  source  with  a time  behavior  of  the  form 
exp(-iut) . The  unknown  $ satisfies 

(5)  4>(o  + s - 4lr  * <r*' 

Equation  (5)  multiplied  by  G minus  equation  (4)  multiplied  by  $ may 
be  integrated  (in  rQ  coordinates)  over  the  domain  of  interest  D 
and  the  resulting  integral  by  Green's  Theorem  reduced  to  an  integral 
over  the  surface  B which  forms  the  boundary  of  D . Tr  any  portion  of 
D is  Infinite  such  as  to  require  a boundary  at  infinity,  appropriate 
limiting  conditions  must  be  Included  for  this  portion  of  B (this  is 
not  difficult  for  most  wave  problems  but  can  require  case  in  some 
cases,  particularly  in  some  water  wave  applications.  This  results  in 
the  Inhomogeneous  Helmholtz  Integral  equation,  e.g.,  Morse  and  Feshbach, 
sec.  7.2  (1973) 

(«  5 GcrWVo 

4 (1/<k)  J [gw  7.sw]*dK. 


8 


i 


with  dA  and  dV  a directed  area  element  and  a volume  element 
o o 

respectively.  The  first  integral  represents  the  contribution  of  the 
inhomogeneous  or  forcing  function  given  in  the  domain  D ; the  second 
represents  the  contribution  of  the  boundary  values  with  the  area  element 
dAQ  positive  away  (outward)  from  the  domain  D . The  gradient  opera- 
tor V is  expressed  in  r coordinates  which  will  be  referred  to  as 
o o 

source  or  integration  point  coordinates;  $(r)  is  evaluated  at  r 
which  will  be  called  the  field  point.  The  parameter  e is  zero  if  the 
field  point  is  exterior  to  the  domain  D (i.e.,  the  delta  function 
6(r  - rQ)  is  zero  everywhere  within  D)  and  is  one  if  the  field  point 
lies  within  the  domain  D . Of  particular  intereet  to  boundary 
Integral  equation  methods  is  the  case  where  the  field  point  lies  on  the 
boundary  B.  For  this  case,  the  boundary  integral  has  a singular  contri- 
bution at  the  point  R » 0 which  may  be  evaluated  by  distorting  the 
actual  boundary  to  exclude  the  field  point  from  the  domain  of  integration 
and  taking  a limiting  process  as  the  distorted  boundary  returns  to  the 
original  form.  For  a plane  boundary  surface,  this  results  in  a value 
of  1/2  for  e ; for  a field  point  at  a corner,  the  appropriate  choice 
for  e is  (l/4ir)  (solid  angle  of  the  domain  excluded),  e.g.,  Banaugh 
and  Goldsmith  (19 63) -Appendix  C.  Furthermore,  the  boundary  integral 
must  be  taken  as  aprincipal  value  Integral,  i.e.,  the  singular  point 
R - 0 is  to  be  excluded  from  the  integration  (since  its  effect  has 
already  been  Included  in  the  modification  of  e to  some  value  other 
than  zero  or  one) . Equation  (6)  may  be  written  as 

(?)  d4>(?)=  eoptikRl/R}  aVo 


•y  s • , • ^ * • . * v • * * 
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where  n is  the  outward  normal  (from  D) . If  the  time  harmonic 
o 

dependence  is  reintroduced,  this  equation  becomes 

+ (i/^Tt } { ,-w/tR  * l/R1] 


B 


These  terms  may  be  thought  of  as  representing  the  inversion  of 
a Fourier  transform  in  time  of  the  time  dependent  functions  iji(r,t) 
and  F(r,t)  which  satisfy: 


(9) 


(f,t)  - (i/c1')  i-vw/n'* 


or  upon  transformation 


(10) 


- 4W  F (?;<->) 


with  tji(r  ;u>)  identified  with  $ and  F(r;o>)  with  f , i.e., 
V\>  (t; u ^ = S it  - 


-CO 


Then 


(11)  £»V(f,t)=  S (^) 


•C.4A*.-*fr-,  1 1 I, 'i 


These  integrals  nay  be  readily  Interpreted  as  tine  retarded  values, 
i.e.,  by  means  of  the  shifting  theorem  in  Fourier  transform  theory. 
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f,  $ and  3$/3nQ  do  not  contain  u>  explicitly  and  the  single  appear- 
ance of  u>  within  the  integral  may  be  interpreted  as  the  transform  of 
a time  derivative.  Thus,  using  a retarded  time  tQ  ■ t - R/c  , this 
Inversion  leads  to 

(12) 

, (1/w)  j ^i(RpT(f„u/^o+  [(4 

which  is  Klrchhoff's  'retarded  time'  integral  equation,  equivalent  to 
the  original  wave  equation.  The  concept  of  a retarded  time  arises 
naturally  in  the  sense  that  an  event  occurring  at  a source  point  tq 
at  time  tQ  is  felt  at  a field  point  r only  at  time  t where 
t - tQ+R/c.  As  mentioned  above,  the  restrictions  on  this  derivation, 
e.g.,  that  <K r,t)  have  a Fourier  transform  zero  initial  values,  limits 
of  integration  over  B which  do  not  depend  on  t etc.,  may  be  removed 
by  a direct  derivation,  e.g.,  Morse  and  Feshbach  (1953)  using  singular 
solutions  or  Friedman  and  Shaw  (1962)  using  a four-dimensional  Green's 
Theorem. 

These  two  Integral  equations,  equations  (6)  and  (12)  form  the  basis 
of  the  boundary  Integral  equation  method.  The  three-dimensional 
Helmholtz  integral  equation  may  be  reduced  to  a two-dimensional  form 
called  the  Weber  integral  equation  without  difficulty; 

\iM  J wTivm  **•>  **• 


(13) 
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from  the  field  point  so  es  to  have  negative  retarded  time  values, 
i.e.,  are  evaluated  at  times  when  the  field  there  has  sero  values  and 
thus  do  not  contribute,  at  least  for  some  Interval  of  time.  This  in 
turn  implies  that  the  limits  of  integration  in  the  integral  depend  on 
time,  increasing  as  time  increases,  i.e.,  a Volterra  Integral  equa- 
tion. Upon  suitable  numerical  approximation,  these  equations  lead  to 
successive,  rather  than  simultaneous,  algebraic  equation,  thereby 
considerably  simplifying  the  solution  procedure.  For  finite  bodies, 
there  will  come  a time  when  all  surface  values  are  nonsero,  but  the 
retardation  still  leads  to  successive  equations.  This  also  presents 
a new  difficulty  however.  The  derivation  of  the  original  Kirchhoff 
equation  was  valid  only  for  continuous  dependent  variables.  A discon- 
tinuous problem,  e.g.,  scattering  of  a step  wave,  requires  additional 
terms.  Fortunately,  these  can  be  evaluated  in  closed  form,  e.g.,  a 
general  discussion  in  Shaw  (1968a)  with  examples  in  Friedman  end  Shaw, 
(1962)  and  Shaw  and  Friedman  (1962). 

An  interesting  sidelight  developed  in  these  last  two  references 
is  the  formulation  of  transient  problems  in  terms  of  a total  rather 
then  a scattered  field.  The  Kirchhoff  retarded  time  equation  as  given 
above  when  applied  to  an  exterior  scattering  problem  must  be  written 
in  terms  of  a dependent  variable  which  decays  appropriately  at  large 
distances  from  the  scatterer  to  allow  a closing  surface  'at  infinity' 
in  Green's  Theorem.  This  clearly  is  the  scattered  (total  minus 
incident)  field.  The  required  boundary  conditions  could  then  be  written 
on  the  scattered  field.  However,  if  such  boundary  conditions  were 
originally  homogeneous  for  the  total  field,  they  would  now  be  inhomo- 
geneous for  the  scattered  field.  This  difficulty  can  be  avoided  by 
writing  an  integral  equation  on  the  incident  field  aa  well.  This  may 
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be  developed  In  several  ways.  The  most  direct  Is  to  consider  the 
boundary  Integral  equation  applied  to  the  scattered  field  at  a boundary 
point.  The  time  harmonic  case  Is  used  as  an  example;  the  approach  Is 
equally  valid  for  the  transient  problem.  The  inhomogeneity,  f(r)  , 

Is  dropped  for  convenience. 


(14) 


w ■ <>«  htr-s*’- wiwfl*' 


Next,  consider  a formulation  for  the  incident  wave,  $w(r)  . This  is 

the  field  which  would  exist  in  the  absence  of  the  scattering  body 

defined  by  the  surface  B.  An  exterior  boundary  integral  equation  on 

* , is  not  available  since  <{>  does  not  in  general  vanish  at  large 
w w 

distances  from  the  scatterer  and  thus  the  integral  over  a closing  sur- 
face 'at  infinity'  in  Green's  theorem  may  not  vanish.  However,  an 
interior  formulation  is  available;  <(>w(r)  “ay  be  expressed  at  some 
interior  point  to  B,  i.e.,  inside  of  the  scattering  p»-face,  by  con- 
sidering the  surface  to  be  absent  but  by  forming  the  integral  equation 
over  those  values  of  $w  and  d^/9n  which  exist  at  those  points 
where  the  scattering  surface  was.  This  leads  to  a similar  boundary 
integral  equation  as  on  4 as  the  interior  point  approaches  the 
location  of  B,  but  with  a change  of  sign  for  the  normal  direction. 
Keeping  n^  the  same  as'  that  in  equation  (9)  leads  to 


{ { 

j 1 $/*•  K 
& 


The  total  field,  $T  ” + ’ may  be  found  by  subtracting  these 

equations,  yielding 


(16) 
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An  alternative,  and  in  some  ways  more  useful,  derivation  may  be  des- 
cribed in  the  form  of  a 'saltus'  problem.  The  exterior  and  interior 
wave  problems  for  the  same  geometry  are  obviously  closely  related.  If 
the  domain  in  question  is  exterior  to  the  boundary  B,  then  e is  one 
for  a field  point  exterior  to  B and  zero  for  a field  point  interior  to 
B and  vice-versa  for  the  interior  problem,  with  a change  in  sign  for 
the  normal,  nQ  , which  is  always  directed  outward  from  the  domain  in 
question.  The  kernel  functions,  G and  3G/3n  , in  the  integral 
equation  formulation  may  be  interpreted  as  layers  of  monopole  and  dipole 
distributions  of  strength  3$/3n  and  $ respectively,  which  replace 
the  boundary  surface  by  an  equivalent  singular  distribution  in  an 
infinite  fluid  domain  including  the  interior  of  the  body  for  the  exter- 
ior case  and  vice-versa,  i.e.,  a 'saltus'  problem.  These  surface 
singularity  distributions  cause  discontinuities  in  3^/dn  and  t of 
strengths  3$/3n  and  $ » i.e.,  they  carry  the  value  of  t and  3+/3n 
from  their  appropriate  values  on  one  side  of  the  surface  B to  zero 
values  on  the  other  side.  This  property  of  distributions  of  monopoles 
and/or  dipoles  is  well  known  and  discussed  in  many  references,  e.g., 
Morse  and  Feshbach  (1953).  Now  reconsider  the  equation  on  , at 
some  field  point  r in  the  fluid  writing  and  3$g/3nQ  within 

the  integral  as  the  actual  discontinuities  in  these  values,  i.e., 
and  td$g/3n]  where  t]  is  the  value  on  the  'fluid*  side  of 
B minus  the  value  on  the  'empty'  side  of  B (that  side  excluded  from 
the  Green's  theorem  integration),  $g  la  replaced  by  and 

, |3*s/fr»oJ  by  tif] , fa4T/3nol  since  *w  is  continuous 
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everywhere  In  the  fluid,  even  across  B. 


(17)  ^ ) ~ ) = ( Vax\) 
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Now  as  the  field  point  r approaches  B,  the  integrals  contribute  a 
tern  l/2(^(r»  leading  again  to  Gq.  (16),  since  4>T  inside  B is 
zero  identically,  by  this  formulation. 

The  introduction  of  the  concept  of  singularity  distributions 
which  replace  the  original  problem  by  an  equivalent  saltus  problem 
leads  naturally  in  a discussion  of  the  ’indirect'  boundary  integral 
formulations;  the  form  above  is  called  the  direct  approach  since  it 
deals  with  the  actual  dependent  variables  4 and  3<|>/3n  'directly'. 

The  indirect  approach  is  based  on  the  observation  that  a distribution 
of  monopoles  or  dipoles  alone  or  in  fact  any  combination  of  these  will 
be  sufficient  to  solve  these  problems;  however,  only  the  direct  approach 
has  the  strengths  of  these  distributions  tied  directly  to  the  actual 
dependent  variables  of  the  problem.  Excellent  discuss-*1'  these 

alternative  formulations  are  given  by  Chertock  (1971)  and  by  Kleinman 
and  Roach  (1974)  in  terms  of  the  present  formulations;  many  other 
classical  discussions  are  also  available,  e.g..  Lamb  (1932) 

Consider  a monopole  surface  distribution  of  strength  a(rQ)  . The 
field  at  a point  r is  given  by 

(18)  $(?)»(WatO  \ f0"^'  J ^ 
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where  a(rQ)  must  be  3$/3no-3$/3no  where  3$/3 nQ  is  the  outward 
normal  derivative  of  the  exterior  (interior)  field  under  consideration 


Is  the  normal  derivative  (in  the  seme  direction)  of  a 


complementary  field;  that  which  satisfied  the  same  boundary  values  for 


4 for  the  Interior  (exterior)  problem.  For  a problem  with  Dirichlet 


boundary  conditions,  this  forms  a 'boundary'  integral  equation  on  the 


monopole  distribution  strength  o(r  ) in  terms  of  the  prescribed 


There  is  no  singular  contribution  as  r approaches  the  sur 


face  B.  Neumann  boundary  conditions  could  be  treated  by  differential 


ing  the  equation  on  $(r)  in  an  n direction  before  approaching  the 


surface;  thus  leads  to  an  alternate  form 


]JcSA.  - 


(J’(ro') 


Similar  results  could  be  obtained  for  the  dipole  layer 


or  any  combination  of  monopole  and  dipole  distributions.  These 


alternate  formulations  have  some  advantages,  as  is  most  apparent  in 


the  solution  of  exterior  time  harmonic  scattering  at  wave  numbers 


near  to  the  eigenvalues  of  the  corresponding  interior  problem.  The 


Fredholm  alternative  indicates  no  unique  solution  to  the  integral 


equation  at  these  values,  even  though  the  exterior  problem  does  have 


a unique  and  well  behaved  solution  for  all  wavenumbers.  This  is  due 


to  the  close  coupling  between  boundary  Integral  equation  formulations 


of  exterior  and  their  complementary  interior  problems.  This  coupling 


has  presented  substantial  difficulty  in  the  application  of  B1E  methods 


to  tine  harmonic  exterior  scattering  problems  as  described  for  example 
by  Chertock  (1971) , Kleinman  and  Roach  (1974)  which  describe  monopole 
and/or  dipole  combinations  to  circumvent  this  problem  and  has  led  to 
several  'modified'  procedures,  e.g.,  those  of  Schenck  (1968),  which 
uses  additional  equations  for  field  points  within  the  scatterer,  Ursell 
(1973)  which  uses  a modified  wave  source  solution  which  satisfies  a 
dissipative  boundary  condition  on  some  surface  interior  to  the  scatterer. 
Of  course,  the  BIE  method  is  quite  suitable  for  the  determination  of 
these  eigenvalues  for  the  interior  problem,  e.g.,  Tai  and  Shaw  (1974), 
DeMey  (1976),  Hutchinson  (1978);  this  is  discussed  further  in  a later 
section. 

Finally,  an  alternative  formulation  to  that  usual  boundary 
integral  equation  needs  to  be  given.  If  3$/3n  is  specified  on  B for 
equations  (6)  or  (16),  for  example,  the  resulting  integral  equation 
is  one  of  the  second  kind  and  numerical  procedures  may  be  used  without 
too  much  difficulty.  If  however,  $ Is  specified,  i.e.,  a Dirichlet 
problem,  the  result  is  an  integral  equation  of  the  fir,_  ..uia  on  the 
dependent  variable,  34>/3n  . While  this  may  be  solved  numerically, 
e.g.,  Shaw  and  English  (1972)  for  transient  scattering  by  a pressure 
release  sphere,  it  may  be  useful  to  convert  the  original  form  into  an 
equation  of  the  second  kind  before  numerical  approximation.  This  is 
done  by  taking  the  field  point  off  of  the  surface  B,  differentiating 
with  respect  to  the  normal  direction  n and  then  allowing  the  field 
point,  at  which  3$/3n  is  now  defined,  to  approach  B; 

(a)  w 4£&*) s * 5 \ '»*>  l * > 
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and  an  analogous  equation  for  the  transient  case  where  such  an  approach 
was  used  by  Shaw  and  Friedman  (1962) . Further  discussion  on  the  dif- 
ficulty of  numerical  solution  of  Fredholm  integral  equations  of  the 
first  kind  is  given  for  example  in  Delves  and  Walsh  (1974). 


NUMERICAL  SOLUTION 


Due  to  the  large  number  of  examples  available  in  literature,  as 
given  in  part  in  the  references  and  bibliography,  and  to  the  several 
review  articles  available,  it  does  not  seem  appropriate  to  give  details 
of  either  numerical  procedures  or  actual  programs.  In  the  sense  that 
this  is  still  a developing  field,  any  of  the  numerical  techniques  and/ 
or  computer  programs  are  subject  to  significant  improvement.  What  does 
appear  appropriate  here  is  a description  of  the  forms  of  numerical 
solution  techniques  and  those  pitfalls  which  have  been  encountered  and, 
where  this  has  been  found  possible,  circumvented. 

Consider  first  the  class  of  Interior  time-harmonic  problems  for 
a finite  region  bounded  by  B where  some  suitable  boundary  conditions 
are  given.  Here  the  homogeneous  and  Inhomogeneous  problems  . ist  be 
distinguished.  The  former  has  as  its  goal  the  determination  of  the 
eigenvalues  and  eigenmodes  for  that  region  with  those  particular  boun- 
dary conditions,  e.g.,  Tai  and  Shaw  (1974),  DaMey  (1976),  Gpegoire, 
Nedelec  and  Planchard  (1975),  Hutchinaon  (1978).  As  an  ex amp Is,  consi- 
der the  two  dimensional  resonance  of  an  acoustic  doatfin  with  a rigid 
boundary,  i.a.,  a Neumann  boundary  condition.  The  governing  homo- 
geneous boundary  lntagral  equation  is 


(22) 


ft 


<^® 
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In  the  simplest  approach,  <J>  is  approximated  by  a piecewise  constant 
over  N segments  of  B,  preferably  equal  in  size,  and  the  field  point, 
rv  , is  placed  at  the  midpoint  of  each  segment  with  *(r_)  for 
r_  ■ r„  defined  as  the  constant  for  that  segment,  . Then 

o K K 

id 

(23)  I 
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with 


(24) 
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The  principle  value  occurs  only  for  that  coefficient  where  the  interval 
L Includes  the  field  point  K , i.e.,  L equals  K . All  other  coef- 
ficients are  calculated  as  ordinary  integrals.  The  integral  for 
excludes  the  point  r equal  to  r ; an  examination  of  the  form  of 
the  integrand  shows  it  to  be  -klT^  (kR)  3R/3no.  As  a first  approxima- 
tion, consider  to  be  evaluated  by  any  numerical 

trapezoidal  rule,  in  which  the  (curved)  segment  is  replaced  by  a 

(large)  number  of  straight  subsegments,  e.g.,  the  intervals  used  in 
the  trapezoidal  rule.  That  straight  segment  containing  the  field  point, 
i.e.,  R equal  zero,  also  has  3R/3no  equal  to  zero  everywhere. 

Since  the  point  R equal  zero  is  excluded  from  the  Integration  (this 
is  a principle  value  Integral) , the  contribution  of  this  subelement  to 
the  coefficient  is  zero.  Of  course,  for  a simple  curved  geometry, 

the  integral  might  be  evaluated  in  closed  form  as  a principle  value 
integral  directly,  but  this  is  unlikely  for  an  arbitrary  shape.  It  is 
useful  from  a numerical  view  to  take  the  subsegments  in  the  integration 
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over  such  thac  Che  actual  field  point  falls  at  the  midpoint  of 

one  of  these  subsegments.  Details  for  the  calculation  of  these 
'influence  coefficients'  (which  are  also  used  in  forced  problems)  can 
be  found  in  almost  any  of  the  references,  e.g.,  Banaugh  and  Goldsmith 
(1963),  for  a somewhat  more  accurate  evaluation. 

The  resulting  system  of  linear  algebraic  equations  must  have  a 
zero  determinant  for  the  coefficient  matrix,  i.e., 

| (VO  ^ 

where  are  functions  of  the  parameter  k , to  have  a nontrivial 

solution.  Those  values  of  k satisfying  this  equation  are  then  the 
eigenvalues  of  the  original  problem. 

Some  care  must  be  taken  with  the  form  of  the  Green's  function 
used.  Some  authors,  deMey  (1977)  and  Hutchinson  (1978),  have  suggested 
that  a simplification  could  be  achieved  by  using  YQ(kR)  in  place 
of  HQ^(kR)  as  the  fundamental  solution.  This  possesses  the  appro- 
priate singular  behavior  and  satisfies  the  Helmholtz  equation  and 
furthermore  leads  to  real  equations  (real  a^)  on  the  (real)  eigenvalues, 

k . Unfortunately,  this  apporach  is  incorrect  since  the  Y (kR) 

o 

fundamental  solution  contains  both  outgoing  and  incoming  waves  whereas 
the  HQ^(kR)  alone  possesses  the  appropriate  behavior  of  outgoing 
waves  at  large  R . This  can  be  readily  seen  in  a simple  example,  a 
circular  boundary,  to  which  the  solutions  are  well  known.  In  this 
particular  case,  described  in  detail  in  the  report  version  of  Tai  and 
Shaw  (1974) , the  boundary  integral  equation  becomes 

us)  (ik)  4,«»  -- 
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where  0 ■ (6-0o)/2  . An  expansion  of  $(0)  as  a Fourier  series 
and  a corresponding  modal  separation  gives 

<^( e->:  £.  [ /w  »iw8t  * 

rrr\-  o 
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These  equations  are  identical;  clearly  either  A or  B could  vanish, 

mm 

but  not  both  for  a nontrivial  solution.  Thus,  in  real  and  imaginary 
parts,  iir  I . 

\ = Jv^djS 

(u  0 = h*  ( T,  UV»w»^  b ' *°-7-  (t  'T" 

o 

The  only  common  root  Is  J'(ka)  equal  to  zero  and  these  are  of  course 

m 

the  correct  colutions  (apart  from  ka  itself  zero) . Mote  however 

that  the  use  of  the  Y (kR)  fundamental  solution  would  have  led  to 

o 

only  the  first  equation  which  has  an  additional  set  of  spurious  roots, 
Yn(ka)  equal  to  zero.  Thus  it  is  necessary  to  keep  the  original  fun- 
, damental  solution  and  pay  the  price  of  the  complex  determinant  calcula- 
tion in  order  to  obtain  the  appropriate  eigenvalues. 

The  corresponding  elgenmodes  can  be  calculated  from  the  original 
set  of  algebraic  equations  by  choosing  a reference  value  for  one  of  the 
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values  4>k  , e.g.,  ^ equal  1.  As  long  as  this  Is  not  a node,  the 

resulting  set  of  N equations  on  N-l  unknowns  may  be  solved  for  the 
eigenmode  on  the  boundary;  Interior  values  may  be  obtained  from  the 
original  Helmholtz  or  Weber  Integral  equation  with  the  field  point 
interior  to  B,  i.e. , 


* 


where  this  is  a simple  quadrature  since  $(rQ)  on  B is  now  known  from 
the  previous  calculation. 

Finally  some  comment  should  be  made  on  the  use  of  boundary  integ- 
ral equation  methods  for  higher  modes.  To  accurately  represent  a rapid- 
ly varying  function,  many  boundary  segments  would  be  required  leading 
to  a large  system  of  algebraic  equations  and  thus  a high  order  deter- 
minant. This  numerical  approach  then  appears  unwieldly;  however  it  is 
for  these  higher  modes  (large  values  of  ka  ) that  the  asymptotic 
solutions  to  the  boundary  integral  equation  may  play  a significant  role. 

If  Dirlchlet  boundary  conditions  are  given,  the  resulting  boundary 
integral  equation  is  one  of  the  first  kind.  These  are  inherently  less 
stable  for  numerical  solution  than  those  of  the  second  kind,  e.g..  Delves 
and  Walsh  (1974) , and  should  probably  be  converted  over  to  equations 
of  the  second  kind,  as  described  above  in  the  section  on  formulation, 
before  proceeding  with  numerical  solution. 

The  interior  inhomogeneous  time  harmonic  problem  has  not  been  as 
widely  studied  as  the  corresponding  exterior  problem,  which  as  great 
application  in  terms  of  wave  scattering  and  radiation  in  our  infinite 
domain.  Indeed,  this  is  probably  the  most  widely  studied  case  for 
the  boundary  integral  equation  method  as  applied  to  wave  problems. 
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The  major  difficulty  with  this  group  of  problems  is  not  the  actual 
numerical  solution  procedure  but  rather  the  inherent  difficulty  in 
obtaining  solutions  to  the  exterior  boundary  integral  equations  at 
wave  numbers  which  match  eigenvalues  for  a corresponding  (but  altered 
boundary  condition)  interior  problem,  e.g.. 

Exterior  Dirichlet  -*•  Interior  Neumann  eigenvalues 
Exterior  Neumann  -►  Interior  Dirichlet  eigenvalues 
This  is  not  a physical  difficulty  inherent  in  the  exterior  problem  since 
there  are  no  'exterior'  eigenvalues  for  these  acoustic  problems,  although 
such  questions  do  arise  for  long  period  water  waves  over  variable  topo- 
graphy, e.g.,  Meyer  (1971).  Rather  this  is  a difficulty  inherent  in 
the  boundary  integral  equation  approach,  as  discussed  by  a number  of 
authors,  with  a corresponding  number  of  methods  of  circumventing  the 
problem,  as  described  in  the  previous  section.  To  illustrate  this  point, 
consider  the  example  of  a two  dimensional  circle  of  radius  a radiating 
into  an  infinite  acoustic  domain  with  a p--“scribed  t-»v.._  harmonic  sur- 
face velocity,  vn  ■ - * vq.  For  simplicity  t-"1"'  the  axisymmetric 

case,  i.e.,  vq  is  constant.  Then  equation  ( t>  ) becomes 

« *,<?)  r (1/0  ) [-AJ.- 
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where  $ is  a constant  $ on  the  surface,  r 

ro 


a , i.e., 
lit 


The  field  point  is  not  placed  on  the  surface  since  the  expansion  of 
the  fundamental  solution,  as  given  in  Morse  and  Feshbach,  sec.  7.2 
(1953),  requires  a different  value  for  |r|  and  |rQ|  . 


v t.x  - r*  cos  (©-©o')! 
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where  e - 1 if  m - 0 and  2 if  m i 0. 
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For  simplicity,  take  0 » 0 and  r - a ; since  tQ  - a,  this  places 
r infinitesimally  less  than  tq  and  also  places  the  field  point 
exterior  to  the  infinite  fluid,  i.e.,  interior  to  the  circular  boundary. 
Then  e ■ 0 and  # 

0=  _ia»s 

, A,  No  \Ao 
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which  agrees  with  the  well  known  separation  of  variables  solution  il 
JQ<ka  ) is  not  zero,  i.e.,  if  k is  not  an  eigenvalue  to  the  corres- 
ponding Dirichlet  problem. 

Apart  from  this  very  real  but  treatable  difficulty,  the  two  and 
three  dimensional  time  harmonic  exterior  problems  are  no  more  diffi- 
cult than  the  interior  counterparts;  this  is  in  strong  contrast:  to  the 
alternative  numerical  methods  of  finite  differences  and  finite  elements. 

Transient  problems  form  the  next  class  for  discussion.  As  mentioned 
above,  the  two  dimensional  transient  fundamental  solution  contains  a 
time  integration  and  it  has  been  found  convenient  to  treat  two  dimen- 
sional problems  in  terms  of  an  extension  into  a three  dimensional  form 


Inal  two  dimensional  body,  e.g.,  Friedman  and  Shaw  (1962).  There  have 

j ' 

been  a number  of  Cwo  and  three  dimensional  problems  considered;  these 
have  been  exclusively  exterior  problems  (see  the  references  and  bibllo- 
graphy) . The  main  advantage  of  time  dependent  boundary  integral  equa- 
tion solutions  is  that  they  lead  to  uncoupled  or  weakly  coupled  systems 
of  algebraic  equations  when  $ is  approximated  by  some  trial  function, 
e.g.,  a piece  wise  constant,  over  the  boundary  elements.  Furthermore, 
solutions  for  discontinuous  forcing  functions,  e.g.,  shock  wave  or 
step  pressure  loading,  can  be  treated  with  the  contribution  from  the 
discontinuity  evaluated  exactly,  e.g.,  Friedman  and  Shaw  (1962),  Shaw 
(1968c) . While  integral  equation  formulations  of  the  second  kind  are 
preferable,  e.g.,  Shaw  and  Friedman  (1962),  some  solutions  have  been 
carried  out  using  integral  equation  formulations  of  the  first  kind, 
e.g.,  Shaw  and  English  (1972). 

As  a simple  illustration,  consider  the  two  dimensional  scattering 
of  a step  plane  pressure  pulse  by  an  infinite  right  a’lgie  wedge  with 
a rigid  boundary  condition;  this  forms  the  'front'  portion  of  the  prob- 
lem described  by  Friedman  and  Shaw  (1962) . The  problem  is  treated 
as  three  dimensional  with  z the  coordinate  axis  In  the  third  direction. 


The  two  surfaces  of  the  wedge  are  taken  as  the  xz  and  yz  planes 
and  the  governing  equation  for  a field  point  at  (x,o,o)  is  then 


, xv  i ( f L <t> cf.Aoi) ciA* 
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Here  S represents  both  surfaces  but  since  3R/3nQ  vanishes  along 
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the  x z surface  (the  singular  contribution  from  the  point  R equal 
o o 

zero  has  already  accounted  for,  i.e.,  the  principle  value  excludes  this 
point)  only  integrations  on  the  yQzo  plane  remain.  Taking  a sym- 
metric incident  wave  for  simplicity,  the  solution  on  the  yo*Q  surface 
is  equal  to  that  on  the  x z surface  and  the  equation  becomes 

i 4>(x(o,ort)s  ^ (*0,0,0  j J { ^ a to  ^dAI° 
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The  limits  of  integration  are  determined  by  the  condition 

(30  **  * 7o 


which  yields 


Xc)  a. 
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Thus,  equation  ( ) has  time  dependent  limits  of  integration.  If  an 

Incident  pressure  field,  pw  » pQ  3<t>/ d t , is  prescribed  as  a unit  step 
function,  H[ct-(x+y)vlT/2] , this  equation  on  the  velocity  potential 
6 may  be  differentiated,  using  Leibnitz1  rule  to  Interchange  the 
spatial  integration  with  the  temporal  differentiation  . . - 4 v* . 

„„  .... 

e » 

The  first  integral  is  evaluated  at  zqu  which  is  Immediately  behind 
the  wave  front;  6 and  36/ 3t  are  known  here  as  the  direct  reflected 
solution;  6 2*  aaro  and  d+/3t  la  2/pq  . Thus  thla  tern  nay  be 
integrated  as  a simple  quadrature,  yielding  a function  A(x,t) . 


The  numerical  solution  of  this  boundary  Integral  equation  may  be 
accomplished  by  taking  p to  have  a constant  value,  p _ , for 

KfL 

(K-l) AY<Yq<KAY  and  for  (L-l) At<t<LAt  . The  time  derivative  may  be 
replaced  by  a backwards  time  difference, 

a-p(o,  . lp(  o,  Vqk,  To, _ Yk. L~1Vl‘l 


The  location  of  the  field  point,  Xj  , enters  through  the  term  R in 
the  retarded  time  t . Define 

* , v ( I iloy)  2*/irt\v  dyodz® 

taw  v,a  = J 1 ,R  ' 

' 0 KW 

k)  if  KIT  r(i/2Tl)J  3C/a/t'o  ol'Ol* 

B KU 

where  Bv  is  the  surface  element  defined  by  (K-1)AY<Y  <KAY  and 
(L-l) At<c<LAt  , such  that  (M-L)At<to<(M-L+l)At  there.  It  is  simpler 
to  use  the  condition  ct  = R in  place  of  equation  (31 ) ; this  would 
have  source  points  contributing  to  the  integral  before  an  appropriate 
time  interval  has  passed  allowing  these  points  to  be  sL^u^-aced  by  tin 
incident  wave  and  then  send  a signal  to  the  field  point.  As  long  as 
the  strength  of  these  contributions  is  zero,  this  integral  in  theory 
does  not  change.  In  practice,  there  is  a small  effect  due  to  the 
finite  discretization  of  the  integral  but  this  effect  is  negligible. 

It  would  not  have  been  appropriate  however  to  make  this  modification 
before  the  wave  front  discontinuity  contribution,  A(x,t)  , was  calcu- 


lated. 


The  boundary  integral  equation  is  then  replaced  by 
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where  is  chosen  at  a 'midpoint',  (J-1/2)AX  and 

(M-l)At  of  the  space  and  time  step  respectively,  with  AY  * AX  ■ ^2  cAt, 
to  allow  the  incident  wave  to  sweep  over  one  space  step  in  one  time  step. 

The  solution  of  this  set  of  linear  algebraic  equations1  is  consider- 
ably simplified  when  they  are  recognized  to  be  successive  rather  than 
simultaneous,  i.e.,  the  value  p^  does  not  appear  on  the  right  hand 
side  of  equation  (36)  since  and  Yj^j  are  both  zero  for  negative 

L and  p , p are  zero  for  zero  and  negative  time,  i.e.,  M <_ L - 1. 

If  t^  were  chosen  to  be  MAt  , i.e.,  the  value  of  p at  the  end 
rather  than  the  beginning  of  the  M*-*1  time  step,  there  would  be  a weak 
coupling  in  that  and  would  both  multiply  values  of  p^ 

on  the  right  hand  side  of  equation  (36).  These  can  be  easily  moved 

over  to  combine  with  the  (1/2)  p term;  the  remainder  of  the  summation 

JM 

does  not  involve  values  of  p at  the  present  time  step.  Choosing 
cAt  < AX  or  AY  should  insure  stability  for  an  explicit  scheme  were 
this  a finite  difference  scheme;  the  appropriate  stability  criteria 
here  do  not  appear  to  have  been  studied.  Clearly  a large  time  step 
(covering  several  space  steps)  would  lead  to  a strongly  coupled  (implicit) 
scheme  while  a small  time  step  (covering  a fraction  of  a space  step) 
would  lead  to  an  uncoupled  system.  The  time  step  and  space  step  used 
in  much  of  this  work  is  matched  to  the  incident  wave  , a step  size 
ratio  of  cAt/AX  equal  to  or  less  than  one  should  be  used  to  remain  within 
the  classical  Courant,  Levy,  Fredricks  stability  criteria,  numerical 
results  are  given  in  Friedman  and  Shaw  (1962)  and  agree  very  well  with 


t^  at  an  'endpoint' 
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an  analytical  solution.  Before  leaving  this  problem  however,  it  is  of 
interest  to  note  that  the  vertex  (o,o)  lies  on  'both'  surfaces  and 
has  dR/dnQ  equal  to  zero  on  both  surfaces,  thereby  causing  the  sur- 
face integral  in  equation  (29)  to  vanish  except  for  the  singular  con- 
tribution at  R ■ 0 which  produces  e equal  to  3/4  instead  of  1/2 
(the  region  excluded  from  the  integration  in  Greens  Theorem  is  3/4  of 
a circle  rather  than  1/2  of  a circle  for  this  boundary  point) . The 
analytical  solution  at  this  vertex  is  then  seen  to  be  $(o,o,t)  ■ 
(4/3)$w(o,o,t)  for  any  form  of  incident  wave.  Furthermore,  boundary 

points  for  r/I  ct  > X > ct  have  a solution  <|>(X,o,t)  ■ 2$  (X,o,t) 

w 

since  the  double  integral  is  zero  for  these  values  even  though  these 
points  have  been  reached  by  the  incident  wave;  this  is  simply  the 
reflection  solution  before  corner  diffraction  effects  reach  the  field 
point. 

The  simple  numerical  cases  described  above  give  the  essence  of 
the  procedures  used.  Improvements  clearly  can  be  made  in  the  form  of 
the  trial  functions  and  boundary  element  shapes  used,  c.fe.,  Cruse  (197-.; 
and,  using  finite  element  methods  to  solve  boundary  integral  equations, 
Jeng  and  Wexler  (1977). 

Alternatively,  Galerkln  or  other  moment  methods  could  be  employed, 
e.g.,  Sharma  (1967).  These  improvements,  though  an  important  part  of 
the  progress  of  this  field,  are  not  in  themselves  necessary  for  a basic 
understanding  of  the  'mechanics'  of  this  numerical  approach. 

APPROXIMATE  SOLUTIONS 

Although  the  emphasis  of  this  chapter,  and  indeed  of  this  text  is 
on  numerical  solution  methods,  this  section  appears  approprlata  in 
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that  it  does  represent  a solution  technique  valid  for  certain  types 
of  problems.  These  problems  generally  fall  in  that  class  which  lends 
itself  to  some  form  of  expansion  in  terms  of  a physical  parameter  which 
is  small.  Typical  of  such  solutions  are  the  Btudles  of  time  harmonic 
wave  radiation  by  shells  forms  of  slightly  nonconcentric  cylindrical 
surfaces,  submerged  in  an  infinite  acoustic  fluid,  Shaw  and  Tai  (1974) 
and  Shaw  (1978a) . Here  an  expansion  of  the  problem  may  be  made  in 
powers  of  the  nonconcentricity  distance;  when  this  distance  is  zero, 
the  outer  and  inner  shell  surfaces  are  concentric  and  the  problem  lends 
itself  to  separation  of  variables.  This  technique  can  be  readily  applied 
to  other  geometries,  e.g.,  spherical,  ellipsoidal,  etc.,  in  which  the 
Helmholtz  equation  separates.  The  details  of  this  procedure  are  left 
to  the  references;  they  basically  involve  an  expansion  of  the  fundamen-r 
tal  solution  into  the  modes  of  the  separable  solution  and  simultaneously 
a further  expansion  into  powers  of  the  separation  distance.  Considera- 
tion of  terms  of  like  order  in  this  (small)  distance  leads  to  a suc- 
cessive system  of  equations  which  are  of  the  same  form  for  each  order 
with  only  a modified  'forcing'  term. 

Another  area  of  approximate  solutions,  which  overlaps  into  the 
use  of  boundary  integral  equation  methods  for  analytical  purposes,  is 
their  use  in  deriving  elastic  plate  or  elastic  shell  theory  from  three 
dimensional  elasticity.  If  one  considers  that  the  ultimate  goal  of 
plate  (or  shell)  theory  is  to  provide  a valid  two  dimensional  represen? 
tat ion  of  an  inherently  three  dimensional  problem,  the  role  of  boundary 
integral  aquation  methods  is  apparent.  They  provide  this  service  without 
approximation.  Since  there  are  in  general  two  distinct  surfaces  (upper 
#nd  lower  or  inner  and  outer)  with  possibly  other  edge  surfaces  over 
which  the  boundary  Integral  is  formed,  it  is  dear  that  the  plate  (shall) 
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thickness  will  enter  Into  the  kernels  (Green's  functions)  through 
their  dependence  on  the  distance  from  field  point  to  Integration 
(source)  point.  An  expansion  of  these  kernels  In  terms  of  this 
thickness  allows  a helrarchy  of  plate  (shell)  theories  to  be  developed 
from  an  essentially  rigorous  foundation  in  which  this  parameter  occurs 
naturally.  Since  this  application  Is  somewhat  removed  from  the  main 
stream  of  this  chapter,  details  are  left  to  a reference,  Shaw  (1978b). 

It  would  be  inappropriate  however  to  leave  this  subject  without  indi- 
cating that  St.  Venants  principle  follows  very  naturally  from  a boundary 
Integral  equation  formulation,  e.g.,  Shaw  (1978c);  indeed,  a remarkable 
proof  of  this  principle  for  the  elastic  half  space  was  given  by  Sternberg 
(1954)  using  a form  of  Betti's  theorem  which  in  turn  is  a basis  for  an 
boundary  Integral  equation  formulation  of  elastostatlcs,  Cruse  and 
Rlsso  (1968,  Cruse  (1968). 

Equally  far  removed  from  this  mainstream  are  the  high  frequency 
approximations  of  optics  and  acoustics,  i.e.,  geometrical  optics 
(acoustics).  Where  these  techniques  are  employed  in  the  treatment  of 
wave  propagation  through  inhomogeneous  media,  boundary  integral  equa- 
tions as  yet  offer  little  competition  in  view  of  the  lack  of  appropriate 
fundamental  solutions. 

However,  boundary  Integral  equation  methods  may  play  a significant 
role  in  the  study  of  high  frequency  scattering  by  complex  structures  — 
especially  in  terms  of  an  expansion  of  the  Integral  equation  kernels 
in  terms  of  the  (small)  parameter,  the  wavenumber.  This  approach 
formed  the  basis  of  much  classical  analysis,  e.g.,  Fresnel-Klrchhoff 
diffraction  theory  (which  contained  other  assumptions  as  well) . 
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CONCLUSIONS  AND  DIRECTIONS 

Probably  the  most  useful  discussion  of  the  application  of  boundary 
Integral  equation  methods  to  wave  propagation  problems  Is  that  of  the 
directions  In  which  these  seem  to  be  heading.  This  Is  best  stated  In 
terms  of  where  these  methods  have  been  used  already  and  conclusions 
as  to  their  present  status.  As  with  most  of  this  chapter,  these  are 
subjective  to  some  extent  but  should  provide  a reasonable  overview. 

In  acoustics,  two  and  three  dimensional  time  harmonic  exterior 
radiation  and  scattering  problems  seem  well  In  hand  e.g.,  reviews  by 
Shaw  (1970),  Chertock  (1971),  Kleinman  and  Roach  (1974).  Primary 
advances  here  would  seem  to  be  in  the  areas  of  Improved  trial  functions, 
since  most  examples  to  date  use  a simple  constant  form,  and  surface 
segmenting  and  in  the  solution  of  the  resulting  system  of  algebraic 
equations  — problems  analogous  somewhat  to  those  currently  found  in 
the  area  of  finite  elements.  It  seems  quite  feasible  that  an  all  pur- 
pose program  can  be  developed  for  arbitrarily  shaped  surfaces.  Some 
care  should  be  taken  however  In  realizing  that  the  strength  of  numeri- 
cal solution  by  boundary  element  discretization  lies  in  the  area 
where  classical  high  frequency  and  low  frequency  procedures  are  inad- 
equate, i.e.,  for  wavelengths  comparable  to  the  length  scale  of  the 
obstacle.  Finally  more  work  Is  required  on  eliminating  those  spurious 
results  which  arise  from  the  complementary  interior  problem  eigenvalues, 
e.g.,  Urcell  (1973),  Schenck  (196B).  In  fact  the  use  of  boundary 
Integral  equation  methods  to  determine  Interior  problem  eigenvalues, 
e.g.,  Tal  and  Shaw  (1974),  DeMey  (1976),  seems  another  area  worth 
further  examination. 

Transient  acoustic  problems  seem  to  also  be  off  to  a good  start 
with  direct  Klrchhoff  retarded  time  integral  equation  Solution  methods 
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suitable  for  early  time,  e.g.,  Friedman  and  Shaw  (1962),  Shaw  and 
Friedman  (1962) , Mltzner  (1967) , and  Fourier  transform  and  corres- 
ponding Helmholtz  Integral  equation  solutions  for  longer  time,  e.g., 
Shippy  (1975).  More  complex  geometries  and  eventually  a general 
numerical  program  for  arbitrarily  shaped  obstacles  seem  to  be  the 
next  steps. 

The  analogous  elastic  problems  have  a further  dilemma.  Should 
their  formulation  be  in  terms  of  displacement  potentials,  which  have 
simple  governing  equations  but  complicated  boundary  relationships, 
e.g.,  Shaw  (1968),  or  directly  in  terms  of  tractions  and  displacements 
which  have  simpler  boundary  representations  but  more  complicated 
governing  equations  (and  thus  kernel  functions),  e.g.,  Cruse  and  Rizzo, 
(1968),  Cruse  (1968). 

Coupled  acoustic  fluid-elastic  structure  problems  seem  a natural 
application  of  boundary  integral  equation  methods,  due  to  their  common 
boundary,  especially  when  the  acoustic  fluid  is  infir^'e  in  extent. 

This  is  true  whether  the  elastic  structure  is  renresenLeu  oy  another 
boundary  integral  equation,  e.g.,  Wilkinson  and  Liu  (1970),  Shaw 
(1973a),  Shaw  (1973b)  or  by  some  other  rperesentation  such  as  finite 
elements. 

Water  wave  problems  form  still  another  class  of  problems  only 
briefly  touched  on  here.  For  long  wave  theory,  the  governing  equation 
is  analogous  to  a two  dimensional  acoustic  wave  equation.  However, 
the  phase  velocity  is  dependent  on  the  water  depth  which  clearly  varies 
in  such  problems  and  thus  the  analogy  is  to  acoustic,  propagation  in 
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an  inhomogeneous  medium.  This  is  further  complicated  by  the  treat- 
ment of  most  of  these  as  exterior  problems,  i.e..  Involving  an  Infinite 
domain.  This  has  lead  to  the  introduction  of  coupled  boundary  Integral 
equation  — other  numerical  methods,  e.g.,  Shaw  (1974c),  9haw  (1975a)jto 
finite  differences,  Shaw  and  Falby  (1977),  Shaw  (1978)  to  finite  ele- 
ments, etc.  for  problems  where  the  infinite  domain  Is  represented  by 
some  constant  depth  (constant  phase  velocity)  at  some  distance  from 
the  local  inhomogeneity.  Indeed,  this  coupling  is  an  area  of  general 
interest,  e.g.,  Zlenklewlcz,  Kelly  and  Bettess  (1977).  For  waves 
over  submerged  obstacles,  a boundary  integral  equation  approach  has 
been  found  useful  based  on  a somewhat  complicated  fundamental  solution 
and  has  been  used  by  Garrison  and  Seetharama  (1971),  Garrison  and  Chow 
(1972),  etc.  with  success.  One  of  the  primary  difficulties  here  is 
that  the  governing  equation  is  simple  — Indeed  it  is  the  Laplace 
equation  — but  the  free  surface  boundary  conditions  are  difficult. 

In  fact,  they  may  be  taken  as  nonlinear  and  still  yield  to  a boundary 
Integral  equation  approach  since  the  governing  differential  equation 
is  linear,  e.g.,  Shaw  (1975). 
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